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Abstract. We review the idea of generating non-extensive stationary distributions based on abstract com- 
position rules of the subsystem energies, in particular the parton cascade method, using a Boltzmann 
equation with relativistic kinematics and modified two-body energy composition rules. The thermodynam- 
ical behavior of such model systems is investigated. As an application hadronic spectra with power-law 
tails are analyzed in the framework of a quark coalescence model. 
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Power-law tailed distributions are abundant in Nature and 
in human technology ranging from high energy particle 
spectra to fluctuations in stock markets or connectivity 
statistics in the Internet. It would be natural to explain 
this abundance by a universal, statistical limiting distri- 
bution since different causes result in similar outcomes. 
A more prestigious attempt is to set such phenomena 
into a united framework of non-extensive thermodynam- 
ics, based on certain generalizations of familiar basic for- 
mulas. In particular generalizations of the Boltzmann - 
Gibbs - Shannon entropy formula were seeked as funding 
stones for such a general treatment [T|[^[5P|^[71[5] . 

Several basic questions arise during this enterprise: 
among those the uniqueness of equilibrium state and the 
entropy function describing irreversibility, the connection 
between composition rules for basic thermodynamical quan- 
tities between two large subsystems and the extensivity 
limit for a system with a large number of degrees of free- 
dom, and the very question that which microscopical mech- 
anisms lead to such a distribution. Is this a sign of non- 
equilibrium, of incomplete equilibrium or just of a new, 
generalized kind of equilibrium? Applying and justifying 
a statistical, least thermodynamical approach to high en- 
ergy heavy ion collisions, as it is being central in the ex- 
perimental quest for quark matter, in particular requires 
clarification of the above questions. Any inference to a 
thermal state and a physical temperature of the quark 
matter from single particle spectra must connect the fit 
parameter measuring the spectral slope to basic princi- 
ples of thermodynamics. 

In recent years we have been succeeding towards an- 
swers to the above problems. After facing the fact that 
transverse momentum spectra fit well to a cut power- 
law distribution towards much higher values than just the 



simple Gibbs-exponential, a particular parton cascade ap- 
proach was suggested in Ref.|9j for generating these dis- 
tributions. It has been observed that the stationary dis- 
tribution generated and maintained by a Boltzmann type 
equation is intimately related to the energy composition 
rule used in two-particle encounters. A simple modifica- 
tion of the kinetic energy addition rule among two part- 
ners, which in high energy collisions is probably related 
to the relativistic kinematics, leads to the observed result. 
A general treatment of abstract composition rules is pre- 
sented in Ref.[10j. where the non-extensivity property is 
related to the deformation of addition rule and hence to 
the deformation of the classical Gibbs exponential. 

The physics' question to begin with is the source of 
non-extensivity, especially for the two most relevant quan- 
tities, energy and entropy. It is relatively easy to con- 
struct examples with non-extensive energy, whenever the 
interaction retails a fractal structure in the phase space 
and therefore cannot be neglected in the large volume - 
large particle number limit, as it is traditional in classical 
thermodynamics. It is much harder to understand non- 
extensive entropy, however. In order to shed some light 
to possible mechanisms by which non-extensivity in one- 
particle variables, like entropy and energy, can occur in 
physical systems, let us investigate a very particular case. 

We assume that in an N-particle system there are two- 
particle correlations left and seek for their relative con- 
tributions to total energy and entropy. For the sake of 
demonstration we regard the following special form of the 
two-particle density: 



Pl2 = fiPl)f{P2)9{ri2), 



(1) 



which is factorizing in the momentum space via one-particle 
distribution functions, but is connected in the coordinate 
space via the pair-distribution function, g{r), of the rela- 
tive coordinates. 
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The trace over states is deterniined via phase space 
integrals, normahzed to satisfy the following conditions in 
d spatial and momentum dimensions: 



Jd'^rl = V, 



(2) 



with V being the total volume, VeS the available volume 
for a partner of a given particle, and n the average (mean) 
density in the system. We normalize the integrals so that 
nV = N and nVcS = N -1. 

Under the above conditions this particular two-particle 
density is normalized to Tr/5i2 = N{N — 1). The partial 
trace over the second particle leads to the familiar one- 
particle distribution function used in kinetic theories: 



Pi=Tr2 (pi2) = (iV-l)/(pi). 

The total entropy of a correlated pair in matter, 
-Tr (pi2 lnpi2) /Tr (P12), is expressed by 



^2 



d^'pi f d-^p2 



{2TTnY J {2TTnY 



d'^n / d'^r2 



P12 lnpi2 
N{N -1)' 



(3) 



(4) 



In calculating this quantity two further individual inte- 
grals occur: 



d'^p 



f{p) ^^f{p), 



ifo 



d'^r g{r) \ng{r) 



Using these notations one arrives at: 



^2 



ifo 



(5) 



(6) 



Generalizing the above expression valid for the two-particle 
density, pi2, to an TV-particle density, P12...N, factorized 
into N{N — l)/2 pair contributions we obtain the follow- 
ing entropy per particle: 



Sn s n 

-rr - = + V^info- 

N n 2 



(7) 



The entropy of such a system is considered to be extensive, 
as long as the specific ratio remains finite in the large 
particle number limit: 



lim — ^ < 00. 



(8) 



In this sense dangerous pair distributions are those, for 
which Vinfo increases with N at fixed mean density n. In 
a familiar piece of matter the pair distribution function 
g{r) approaches one at large distances, in these cases Vinfo 
is finite and hence the entropy is extensive. In case of 
a quark gluon plasma, however, some infrared magnetic 
modes remain non-perturbative and hence long range cor- 
relations remain. As a consequence —glng may not tend 





Fig. 1. Schematic plots of the pair distribution function and 
the corresponding energy and entropy contributions in short 
(top) and long range (bottom) correlated matter assuming a 
1/r and a ar type pair potential, respectively. 



to zero fast enough and therefore the integral Vinfo in 
eq.(l5]) may increase as a function of N. For example con- 
sidering power-law type pair distribution functions, like 
g{r) — r°/(l -|- r^), the corresponding integrals up to a 
large radius, R scale like Ves ~ Jl^+a-b ^^^^^ jjj^g y:^^^ ^ 
Jld+a-b -^j^ In this case the specific entropy for large N 
becomes 



N 



— + const. (a 

n 



b)NlnN 



(9) 



with some unspecified constant. For a ^ b this would lead 

to a non-extensive entropy. For other possible sources of 

non-extensive entropy see Ref . pT] . 

The total energy can be calculated in a similar way. 

Assuming a v{ri2) pair-potential depending on the relative 

coordinate only and individual kinetic energies, K{pi), we 

arrive at _ _ 

En e n. 



N 



with 



d-^p 
{2nhy 



f{p)K{p), 



(10) 



Fpair = d'^r g{r)v{r). 



(11) 



Some typical g{r) functions are shown in Fig[TJ On the 
top figure a pair distribution function tending to one at 
large distances and a Coulomb-like pair potential, while 
on the bottom figure a linear confining potential, v{r) ^ 
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r and a power-law tailed pair distribution function are 
assumed. 

In most physical systems studied traditionally in ther- 
modynamics, like gases, liquids, plasmas, etc. the function 
g{r) approaches the value one at large distances. Therefore 
there are no non-extensive contributions to the entropy 
per particle. For the energy the situation is different so 



far, since v{r) 



might not approach zero for large 



distances fast enough, producing this way a contribution 
to the energy per particle, Em /N which may even diverge 
in the large N limit. Such a case is an unscreened l/r-like 
potential in three dimensions. 



2 Additivity, extensivity and abstract 
composition rules 

Non-extensive quantities, whose amount per particle is not 
finite in the thermodynamic limit, are also not additive, 
because the repeated composition by simple addition rules 
always leads to a result proportional to the number of 
steps. It is possible, however, that one is able to find an- 
other quantity, a certain function of the non-additive one, 
which is additive. This way the non-extensive thermody- 
namics can be treated by mathematical algorithms which 
were designed for additive composition rules. In this sec- 
tion we analyze the mathematical background of compos- 
ing energy and/or entropy of subsystems and then repeat- 
ing this composition. The thermodynamical limit is ap- 
pointed to the infinite repetition of the composition with 
an infinitesimal amount 1101. 



2.1 General rules and thermodynamical limit 

Let us denote an abstract pairwise composition rule by the 
mapping (x,y) — » h{x,y). The important question arises, 
that what happens if we repeat such a composition rule 
arbitrarily long, each time applying to an infinitesimal 
amount: This way one deals with the thermodynamical 
limit of composition rules corresponding to the energy or 
the entropy. The effective rule in this limit, which applies 
to results of repeated rules themselves, has special prop- 
erties then. 

From the starting rule we demand only a trivial prop- 
erty: that the composition with zero should be the identity 



h(x, 0) = X. 



(12) 



We do not assume in general symmetry (commutativity) , 
such as h{y, x) = h{x, y) nor we demand associativity 



h{h{x,y),z) = h{x,h{y,z)). 



(13) 



Her we note that the general solution of the associativity 
equation (jl3p is given by 



h{x,y) ^ X-\X{x) + X{y)) 



(14) 



with X{x) being a strict monotonic function [12] . We shall 
refer to this mapping function as the " formal logarithm" , 



because it maps the arbitrary composition rule /i(x, y) to 
the addition by taking the X-function of cq. lM]) : 



X{h{x,y))^X{x)+X{y). 



(15) 



Due to this construction the generalized analogs to classi- 
cal extensive (and additive) quantities are their formal log- 
arithms, whenever the composition rule is associative. As 
a consequence stationary distributions, in particular those 
obtained by solving generalized Boltzmann equations [9], 
are the Gibbs exponentials of the formal logarithm, 



(16) 



Let us now regard a large number of iterations, TV, of 
a general composition rule. We apply it to a small amount 
yjN and repeat this (N — 1) times, constructing this way 
the quantity 



N-\ 



We consider the large- limit. 



N 



lim XN{y) < oo, 



(17) 



(18) 



if this is finite for a finite y, we can apply all formulas 
of classical thermodynamics usually applied to extensive 
quantities. Such a limiting quantity is extensive, but not 
necessarily additive. Our goal is to obtain the asymptotic 
composition rule. 



(19) 



in the limit iVi, N2 — > 00. The recursion for the n-th step 
of this repetitive composition is given by 



(20) 



starting with xq = 0. Subtracting a;„_i — /i(x„_i,0) from 
both sides we arrive at 

y 

Xn - Xn-1 = h{Xn-l, — ) - h{Xn-l,0). (21) 

Denoting by t = (n — l)/iV the extensivity share already 
achieved, one step takes At — l/N, and the above recur- 
sion can be Taylor-expanded for a small y/N^ yAt: 



d 

x{t+At)~x{t)^yAt —h{x{t),y) 
oy 



-0{At^). (22) 



y=0+ 



In the large N {At 0) limit this becomes equivalent 
to a differential equation similar to a renormalization flow 
equation: 

^^yh',{x,0+). (23) 

In this expression h'2{x,0^) denotes the partial derivative 
of the rule h{x,y) with respect to its second argument 
taken when this value approaches zero from above. Note 
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that the uniformity of subdivisions to y/N is not neces- 
sary; aU infinitesimal divisions summing up to y by t = f 
lead to the same differential flow equation. 
The solution of eq. ([23l) . 



L{x) 



dz 



h',iz,0+) 



yt, 



(24) 



defines the additive mapping of x, i.e. the formal loga- 
rithm L{x). By the help of this the following asymptotic 
composition rule arises: 



Xl2 



'f{xi,x2)^L-^{L{xi)+L{x2)); (25) 



it is already associative and commutative. Commutativity 
is trivial and associativity is also easily proved: 

<^(<^(a;i,a;2),a;3) = L^'^ {L{ip{xi,X2)) + Lix^)) 
^ L-^ {L{xi) + L{x2) + L{X3)) 
^ L-^ {L{xi)+L{ip{x2,X3))) 
= 'p{xi,f{x2,X3))- (26) 

It is interesting to check that all associative rules are 
mapped to themselves in the above limit. Given an asso- 
ciative composition rule, h{x, y), it possesses a formal log- 
arithm, X{x), which is additive: 



X{h{x,y))^X{x) + X{y). 



(27) 



Now taking the derivative of this equality with respect to 
the second argument we obtain 



X'{h)dh/dy = X'{y) 
which taken at y = becomes 

X'iO) 



h'^{x,0^ 



X'{h{x,0)) 



(28) 



(29) 



Due to the property h{x,0) = x (equivalently X{0) — 0) 
the formal logarithm of the asymptotic composition rule 
is given by 



^ ' ' X'{0) X'{Q)' 



(30) 



it is proportional to the formal logarithm of the starting 
rule. Therefore the asymptotic rule is exactly the same as 
we begun with: 1^9(2:, y) = h{x, y). The freedom in a factor 
of the formal logarithm is used to set ^'(0) = 1. This way 
any associative composition rule describes a limiting rule 
of a class of non-associative rules. 



2.2 Deformed logarithms and deformed exponentials 

The stationary distribution eg.pB]) in the large- limit 
contains the formal logarithm, L{x). In fact the composed 
function, ea = exp oL is the one, which is frequently called 



a 'deformed exponential' in the literature. Its inverse, In^ = 
L^^ o In is then the corresponding 'deformed logarithm'. 
These functions are inverse to each other. Further proper- 
ties of the traditional exponential and logarithm functions 
are, however, not automatically inherited. In particular 
reciprocals and negatives follow different rules as we are 
used to. 



In the particular case of scaling formal logarithms, 

1 



La{x) = -Li{ax) 



(31) 



several interesting identities hold, among others the fol- 
lowings: 



x, 
1 



Lo{x) 

liia{l/x) = -ln_a(a;), 
l/ea{x) = e^ai-x) 



(32) 



Since a = q — 1, the a* — —a duality corresponds to the 
q* — 2 ~ q Tsallis-duality. This can be important for the 
particle-hole relation for fermions: 



1 - 



1 



1 



ea{-x) + l e„a(a;) + l' 



(33) 



Let us now list some important particular rules and 
their asymptotic pendants considered in applications of 
non-extensive statistics to physical systems. 

The trivial (and classical) addition is the simplest com- 
position rule: h{x,y) — x + y. In this case /i2(x,0+) — 1 
and one obtains 



L{x) — J dz 



(34) 



The original Gibbs exponentials, e^'^^/Z, result as sta- 
tionary distributions from any Monte Carlo type algo- 
rithm using the additive composition rule. The asymptotic 
rule is also the addition ip{x, y) — x + y. 

Another rule leading to the so-called g-exponential dis- 
tribution |J3j is given by h{x,y) — x + y -\- axy with the 
parameter a proportional to 5 — 1 occurring in the Pareto- 
Tsallis distribution. Now one obtains h2{x,0~^) =^ 1 + ax 
and 



Lix) 



dz 



1 



1 



ln(l 



(35) 



This formal logarithm leads to a stationary distribution 
with power-law tail as the function composition exp o L 
on the power — /?: 



1 



;i + aE) 



-13/a 



(36) 



On the other hand, assuming such a non-additive compo- 
sition rule for the generalized entropy, a special formula 
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can be constructed as the expectation value of the inverse 
of this function, of the deformed logarithm, L^^ o In. One 
obtains 



S 



f- = - {.f-^-f)- 



(37) 



The asymptotic composition rule again coincides with the 
original one: (/^(x, y) = x + y + axy. We note here that 
the formal logarithm of the integrated expression is the 
(additive) Renyi entropy: 



L{S) 



1 



In / P, 



(38) 



with a = 1 — g and J f = 1. 

A further rule has been suggested by KaniadakisfT^, 
based on the sinh function. The formal logarithm is given 
as 

L{x) = -Arsh(Ka;), (39) 

K 

and its inverse becomes L^^{t) — smh(Kt)/ k. The station- 
ary distribution, composed by exp o L, is 



1 



(40) 



For large arguments it gives a power-law in the momen- 
tum p and hence also in the relativistic energy. The corre- 
sponding entropy formula is the average of o In over 
the allowed phase space: 



S 



K 



— sinh(Klnf)= / . 



The composition formula can be reduced to 



h{x, y) 



x\J 1 + K^y^ + yyjl + n^x'^ . 



(41) 



(42) 



For low arguments it is additive, h{x, y) x + y, for high 
ones it is multiplicative, h(x,y) « 2Kxy. It has been mo- 
tivated by the relativistic kinematics of massive particles. 
Interpreting the parameter as k = l/mc, one deals with 
Kp = sinhry, so the formal logarithm becomes proportional 
to the rapidity, L{p) — mcq. This implies a stationary dis- 
tribution like exp{—pmcr]), which has not yet ever been 
observed in particle spectra stemming from relativistic 
heavy ion collisions. For such a purpose it is tempting to 
consider some further scenarios based on other quantities 
than suggested above (see next section). 

The rule leading to a stretched exponential station- 
ary distribution, often considered in problems related to 
anomalous diffusion and Levy-flights, is given by h{x, y) = 

{x^ + y^y^^^ . Here the partial derivative is evaluated at a 
small positive argument, e = y/2N. One obtains /i2(a;, e) = 
c{e)x^~^ with a factor depending on e and for given val- 
ues of h diverging in the e = limit. However, this can 
be accommodated by our procedure; we obtain the formal 
logarithm L{x) — c[e)x^ /b, and therefore the asymptotic 

rule ^p{x, y) — {x^ + y'')^^^. Again, constant factors in the 



formal logarithm can be eliminated without loss of any 
information. 

Now let us investigate a non-associative rule; its asymp- 
totic limit cannot be itself. We regard a linear combination 
of arithmetic and harmonic means: 

xy 



h{x, y) = X + y + a- 



V 



(43) 



The rescaling flow derivative is given by /i2(a;, 0+) = 1 + a 
and ~ being a constant - it leads to L(x) = x/ (\ + a) and 
with that to the addition as the asymptotic rule: f{x, y) — 
x + y. 

As an interesting rule we discuss the relativistic for- 
mula for collinear velocity composition. 



h{x,y) = 



y 



1 -I- xy/c^ 



(44) 



This rule is associative, and it also preserves its form in 
the thermodynamic limit. The fiducial derivative is given 
by h2{x, 0+) = l — x^/c^ and the formal logarithm, L{x) = 
catanh (x/c) turns out to be the rapidity. The asymptotic 
composition rule recovers the original one. 

There are also general types of composition rules, which 
mutate into a simpler asymptotic form. For our discussion 
particularly important are rules of the form 



h{x,y) = x + y + G{xy) 



(45) 



with a general function G(z), restricted by the property 
G(0) = only. In this case h'2{x,Q) = 1-1- G'[Q)x asymp- 
totically leads to a Tsallis-Pareto distribution with the 
parameter 5 — 1 = G"(0). 

2.3 Extreme relativistic kinematics 

In this section we review a particular type of pair inter- 
action, which can be expressed as a function of the ki- 
netic energies of the individual particles. The relation to 
relativistic kinematics is established by the fact, that we 
consider such dependence through the Lorentz-invariant 
relative four-momentum square variable: 



E12 = Ei+E2 + UiQ^). 



(46) 



We study whether relativistic speeds alone can cause "non- 
extensivity" , i.e. a power-law tailed kinetic energy distri- 
bution. The relativistic formula for is given by: 



{Pi - P2f - [El - E2f 



(47) 



with Pi,Ei being relativistic momenta and full energies 
of interacting bodies. Expressed by the energies and the 
angle between the two momenta this becomes a linear 
expression of cos 0: 

^2{EiE2~PiP2Cos0)-{ml + ml) (48) 



with Pi ~ ^ El — nrf for i = 1,2. Here we use relativistic 
units (c — 1) and assume the masses mi and m2, respec- 
tively, for the interacting partners. It is useful to note that 
writing eq. (|48)) as = 2{A — Bcos0) we have 



A±B = E1E2 



1 



±PlP2- 



(49) 
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For the sake of simphfication we average over the rel- 
ative directions of the respective momenta and obtain 

(C/(Q2)) = ^ J U{2A-2Bcos0)sm0de 



F{2A + 2B) -F{2A-2B) 
AB 



(50) 



with U (w) = dF/dw. It is easy to derive by the substitu- 
tion w = 2{A — B COS0). The rule for the kinetic energy, 
Ki = Ei — rrii , composition is given by 

F(2A + 2B) - F(2A-2B) , , 
Kx2 = Ki + K2 + '—^ (51) 

The quantities A and B^ can be expressed by the respec- 
tive kinetic energies and masses: 

1 



A = K1K2 + {m2Ki + miK2) - ^(™i - "^2) 
52 = KiK2{Ki + 2mi){K2 + 2m2). 



(52) 



One observes that the product of kinetic energies occurs 
due to kinematic reasons. 

Taylor expanding the integral of the unknown function 
U{w) around w = 2A and ensuring the h{x, 0) — X, as 
well as the h{0,y) — y property, we obtain the following 
composition rule for the relativistic kinetic energies: 



h{x, y) ^ x + y ~U {2m2X + 11112) - U (2miy + TO12) 

(4^2 )J 



+C/(mi2)+;^{/(2jn2A) 



(2j + 1)! 



(53) 



with mi2 — — (toi— m2)2, A — xy + {m2X + miy) + TO12/2 
and AB^ — Axy{x + 2mi){y + 2m2)- For unequal masses, 
TOi ^ TO2 this composition rule is not symmetric. Since 



dB' 
dy 



dA, , 

— (a;,0) = mi+x, 
dy 

(x, 0) — 2m2x{x + 2mi), 



(54) 



the derivative leading to the formal logarithm of the asymp- 
totic rule becomes an expression with a finite number of 
terms 



h'2{x,0) = l-2mi U'{mi2) + 2{mi+x) U' {z) 
4 

+ -m2x{2mi+x)U"{z), 



(55) 



with z = 2A{x,0) — 2m2X + TO12. In all traditional ap- 
proaches the interaction energy U is independent of Q^. 
In such cases h'2{x,0) — 1 and the simple addition is the 
asymptotic composition rule. Therefore the stationary en- 
ergy distribution is of Boltzmann-Gibbs type. For de- 
pendent interactions on the other hand it is important to 
consider the extreme relativistic kinematics. In this case 
the replacement mi = m2 = leads directly to 



h2{x,0) = l + 2xU'{0). 



(56) 



As discussed in the previous subsection this generates a 
Tsallis-Pareto distribution in the relativistic kinetic en- 
ergy. This result includes for U' = the traditional mo- 
mentum independent interaction case leading to the ad- 
dition as asymptotic rule for non-relativistic kinetic en- 
ergies, and hence to the Boltzmann-Gibbs distribution. 
We note that in the relativistic kinematics the linear as- 
sumption, U' = a = const, also leads to a Tsallis-Pareto 
distribution due to h'2{x, 0) = 1 + 2ax. 

2.4 Generalized entropies to each composition rule 

There are two possible approaches in constructing a gener- 
alized entropy formula: i) either to use a non-additive en- 
tropy for independent events with factorizing probability, 
or ii) to search for an additive entropy while the common 
probability is not factorizing in the individual probabili- 
ties. In both cases the entropy density function, a{p) to a 
probability p can be obtained from the composition rule 
h{x,y). 

First we consider a non-additive entropy formula for 
factorizing probabilities, i.e. 

w,j (T{w,j) = hyY p, cr {p, ) , ^ cr (9j ) (57) 

with Wij = PiQj. We would like to construct the function 
cr(p) by knowing h(x,y). Let us inspect the equipartition 
case, Pi = l/Ni, qj = 1/N2- In this case Wij — l/{NiN2). 
Eq.dSIl leads to 



a{ab) — h ((7(a), <j{b)) 



(58) 



with a — 1/Ni and b — 1/N2- This requires the same 
composition rule for micro- and macro-entropy: 



Ypiqjh{(T{pi),a{qj)) = h ^pjcr(pj), ^ 9jCr(gj) 

ij \ i 3 

(59) 

This /i-extensivity can so far only be satisfied by the Tsal- 
lis rule h{x^ y) = x -\-y -\- axy. On the other hand if go = 1 
and all other qj — for j ^ 0, we obtain two constraints: 



ct(p») = h{a{pi), (t{1)) 
a{Q)^h{a{pi),n{Q)) 



(60) 



from which it follows h{x, 0) = x with a{l) = (the un- 
expectedness of a sure event is zero) and cr(0) — 00, too. 

Based on the properties of the known h{x,y), in the 
thermodynamical limit it is associative and hence possess 
a formal logarithm, L{x). Therefore 



L{a{ab)) = L{a{a))+L{a{b)), 



(61) 



whose general solution is given by i o a = /3 In. According 
to the tradition /3 = — 1 in units of the Boltzmann con- 
stant, As = 1, and therefore the entropy density function 
is expressed by the deformed logarithm: 



aip)^L-'i-\np)=lna - 



(62) 
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It is possible to ask another question: if the construc- 
tion rule for the common probability is not the simple 
product, but it is known, what should the entropy density 
function be in order to lead to the addition rule for the 
total entropy. So given the formula 



w,i = e 



h(\n Pi ,ln Qj ) 



(63) 



how to construct a{p) such that 



is fulfilled. We have 



c = Wii — e 



h(\n a, In b) 



(65) 



as the known composite probability and seek for the en- 
tropy density function, cr(a) satisfying 



ca{c) = ab {cr{a) + a{b)) . 



(66) 



We solve this functional equation by deriving with respect 
to b and take the result at & = 1. Since 



|=.e''('->'')/.',(lna,ln&)i 



(67) 



we arrive at 



a/i2(ln a, 0) (cr(a) + a(T'{a)) — aa{a) + aa{l) + aa'{l). 

(68) 

Using now that cr(l) — and h'2{x, 0) — 1/L'{x) with the 
formal logarithm L associated to the composition rule, h, 
we obtain - using the variable x = In a - 



ax 



(69) 



with P = (t'(1) constant. The final solution is expressed 
by the formal logarithm of the asymptotic rule as 

/•In a 

cr(a) = /3e-^('" / L'(u)e"-^(")du. (70) 



It is interesting to note, that using L^^ of the deformed 
logarithm as the function L belonging to the product com- 
position rule (|65|) . one assumes hia{wij) — Ina(pi) -I- hia{qj) 
and arrives at 

pa(p) = -e'""P r'^\-^ea{u)du. (71) 
Jo 



3 Non-extensive Boltzmann equation 



It is a false belief that only the exponential distribution 
can be the stationary solution to the Boltzmann equa- 
tion; this statement is true only i) if the two-particle dis- 
tributions factorize, ii) the two-particle energies are addi- 
tively composed from the single-particle energies (-E12 = 



h{Ei, E2) — El + E2 ) and iii) the collision rate is multi- 
linear in the two-particle (and two-hole) densities. A gen- 
eralization of the original Boltzmann equation has been 
pioneered by Kaniadakis[TS] investigating nonlinear den- 
sity dependence of the collision rates. An "H'^ theorem for 
the particular Tsallis form of the collision rate has been 
derived by Lima, Silva and Plastino[16|. 

A possible generalization of the Boltzmann equation 
uses an altered form of the 'Stosszahlansatz' and allows for 
an evolution equation of a function of the original phase 
space occupation factor, F{f): 



DF{h) = j W1234 (G34 - G12) 



with 



234 



DF 



P 



(72) 



(73) 



total (Vlasov-) derivative, with a 1234-symmetric colli- 
sion rate including Dirac-delta distributions for momen- 
tum and energy composition rules in two-to-two collisions 
(which also may be of generalized type by using corre- 
sponding formal logarithms), and finally the generalized 
product for the two-particle density factor, 



G12 = 6, (ln,(/i)+ln„.(/2)) 



(74) 



using the deformed exponential and logarithm functions. 
Based on this, a particular expression for the entropy cur- 
rent density can be defined: 



(75) 



The entropy density form, <7{F) always can be constructed 
in a way, that the second theorem of thermodynamics is 
fulfilled. The local source for the entropy is namely given 

by 

d^S^ I a\F{h))DF{h). (76) 



Utilizing the generalized Boltzmann equation (|72p and ex- 
changing the index 1 with 2, 3 and 4 while z/;i234 stays 
invariant and obviously Gy = Gji, one arrives at 

9^5''' = \ j "^1234 (ct'i +(^'2-^3- 0-4) (G12 - G34) 

1234 

(77) 

with cr- = cr'(F(/,)) for i = 1,2,3,4. This quantity is 
always non-negative, i.e. 

(<?(Gi2) - <Z'(G34)) (G12 - G34) > 0, (78) 

if and only if 

<Z'(Gi2) = cj\F{h)) + a'(F(/2)) (79) 

is a monotonic rising function. Inspecting the generalized 
Stosszahlansatz eq. (f74| one finds that this splitting to the 
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sum of respective functions of /i and /2 is only possible, 
if (^{t) oc Ina(i). Therefore we conclude that 

a'{F{f))=a\nM)+P (80) 

with a > and /3 undetermined constants. (This deriva- 
tion followed the spirit of Ref.|15j.l 

The generalized entropy density as a function of the 
one-particle phase space occupation density is hence given 

by 

^if)= J F'{f){a\na{f) + P)df. (81) 

The traditional Boltzmann formula arises for F{f) = f 
and lna(/) = lii(/) (i.e. a — 0). Lavagno et.al. [17] consid- 
ered F{f) = and ln,(/) = (/«-i - l)/{q - 1) (i.e. a = 
{q — 1) and Tsallis composition rule for In/). In the case 
of h{x, y) = x + y + axy one considers lna(a;) = (a:° — l)/a, 
ea{t) = (l + at)^/" and Gu = (/f + /2 - 1)^/°- For a small 
a parameter it is G12 ~ /i/2(l — aln(/i) ln(/2) + ...). 

We note that the detailed balance distribution is given 
by the condition G12 = while the corresponding en- 
ergy composition rule apphes L{Ei) + L{E2) = L{E'i,) + 
L{Ei). This is possible only if \na{fi) = -{L{Ei) - ^i)/T, 
so 

r{E) . e. (^^^) . (82) 

The parameters T and /i are arbitrary constants for being 
a stationary solution of the generalized Boltzmann equa- 
tion, but they can be related to the total energy and par- 
ticle number in a given application. 



3.1 Deformed energy composition rules in parton 
cascade 

Our fist numerical approach [9] was restricted to the use 
of abstract composition rules in the energy balance part: 
we equated the energy of the reacting parts before and 
after the collision via an abstract energy composition rule 

h{Ei,E2) = h{E3,E^). (83) 

Although this rule cannot be specified without further 
knowledge, according to our results presented in the previ- 
ous section, in the thermodynamical limit an asymptotic 
rule can be considered, with a formal logarithm. The par- 
ton cascade simulation based on a Boltzmann equation is 
hence modified by considering 

L{Ei) + L{E2) = L{E:i) + L{Ei). (84) 

At the same time we applied F{f) — f and a = 0. Ap- 
plying such a general energy composition rule considered 
in the thermodynamical limit, the rate of change of the 
one-particle distribution is given by 

A = / W1234 [fsU - /1/2] . (85) 
^234 



Tsallis kinematics E = 4 GeV, P = 3.6 GeV, m = 0.14 GeV 




a = +0.0 

a = +0.1 

a = -0.1 

a = -0.2 



-2.5 -2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 
Tsallis kinematics E = 4 GeV, P = 2 GeV, m = 0.94 GeV 







/ // / V 




A \ \ P3\^ 









-2 -1.5 -1 -0.5 0.5 1 1.5 2 
Px 

Fig. 2. Momentum vectors for pions (upper) and protons 
(lower) a.t E — 4 GeV pair energy according to ea. (|86|l with 
the rule h{x,y) = x + y -\- axy for the energies. 

with the symmetric transition probability W1234 including 
the constraint 

A = 5\pi+p2-P3-Pi)5{h{EuE2) ~ hiE3,Ei)). (86) 

In the figure [2] the possible pairs of momentum vectors 
are shown for the h{x,y) = x + y + axy energy compo- 
sition rule for pions (m = 0.14 GeV mass) and protons 
(m = 0.94 GeV mass), respectively. The two-dimensional 
cuts for the endpoints of the respective vectors form an 
ellipsoid in the traditional a — Q case, while this surface 
is deformed for nonzero extensivity parameters, as seen in 
the figure. 

In a stationary state the f{Ei) distributions depend on 
the phase space points through the energy variables only 
(this is to be checked on experimentally observed hadron 
transverse momentum spectra at mid-rapidity by the so 
called TTiT-scaling) and the detailed balance principle re- 
quires 

f{E,)f{E2) = f{E3)f{E,). (87) 

With the generalized constraint (|86p this relation is satis- 
fied by 

f{E) - f{0)exp{^L{E)/T). (88) 

For the Tsallis- type energy addition rule[l8lfT9|. one ob- 
tains cut power-law stationary distribution, 

fiE) = f{0)il + bEr'^'^ . (89) 
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Connecting this to the Tsalhs paranietrization one uses 
q = 1 + bT. Since the energy addition rule conserves in a 
two by two collision the quantity h{Ei, E2), the new en- 
ergies after such an event lie on the /i(£'3, i?4)= const ant 
surface. Due to the additivity of the formal logarithm of 
the single particle kinetic energies, L{Ei), the total sum 
Ltot = '^iL{Ei), is a conserved quantity. This rule was 
applied in numerical simulations [^122] . During the numer- 
ical searches for stationary distributions only the tacit as- 
sumption of constant transition probability rates has been 
applied; the evolution results are obtained in terms of the 
number of pairwise momentum exchange events, not in 
terms of real time. 



3.2 Random momenta 

Parton cascade simulations usually consider pairwise colli- 
sions with energy and momentum conservation inside the 
two-particle system. The pairs to collide are chosen ran- 
domly from an ensemble of particles and the new mo- 
menta are generated randomly according to the above 
constraints. This way the probability is uniform in the 
two-particle phase space, provided the conditions for mo- 
mentum and energy sums (in our more general case for 
the energy composition) are satisfied: 

d?w = wo5{pi +P2- P)S{h{Ei,E2) - H)d'vxd^Vi- (90) 

The constant wq is fixed by the normalization of the in- 
tegral of this probability density to one (or to the actual 
collision rate in real-time simulations). Since there are six 
degrees of freedom and four constraints, two free quan- 
tities have to be chosen randomly. It is, however, a del- 
icate procedure to ensure the random uniformity in the 
two-particle phase space for a general energy composition 
rule. 

It is customary to introduce the sum and difference of 
the momentum vectors by 



Pi,2 = \^P±q. 



(91) 



Using this notation the momentum sum constraint can be 
integrated out trivially and - since the Jacobean of the 
transformation (j9ip is one - we arrive at 



wo5{h{Ei,E2)-H)d^q. 



For the addition rule, h{Ei,E2) = Ei + E2, it is enough 
to obtain the direction of the vector q accordingly while 
its magnitude is constrained by the energy sum. It is a 
straightforward task to do it in the center of mass sys- 
tem, where the momentum sum vector, P, vanishes: the 
direction of the difference vector q in this system is uni- 
form on a spherical surface. A Lorentz-transformation into 
this system, a random azimuthal angle and a random co- 
sine, and finally a back transformation provide the new 
momenta after a collision. 

Since we are dealing with a constraint more general 
in the energy variables, first we transform the problem of 



randomly choosing the difference vector q into a problem 
of choosing proper energies after the collision. The energies 
are expressed by the free dispersion relations 



El2-ml 2 - ^P^ + q^±Pqcost 



(93) 



where P and q denote the lengths of the corresponding 
vectors and 9 the angle between them. From this two equa- 
tions one easily derives the following energy differentials: 

2EidEi = 2qdq + P cos Odq - Pq sin 9de, 
2E2dE2 = 2qdq - P cos 0dq + Pq sin OdO. (94) 

The phase space volume element can be expressed easily 
by using the wedge product form: 

d^q = dqAq sin OdO A qdc/) (95) 

which upon using eq. ()94p can be written as 



,3 E1E2 
d 9= — ^ 



dEi AdE2A> 



(96) 



Now using the energy composition constraint we arrive at 
a probability density which is not uniform in the energy: 

d^w = WQ6{h{Ei,E2) - H)^^dEidE2d(p. (97) 

One uses the constraint to eliminate say E2 from the above 
formula and considers 



d^w ~ Wo 



E1E2 



PK{Ei,E2 



-dEic 



(98) 



In the general case the differential probability density, 
dw/dEi, is a complicated function of the energy. Its inte- 
gral, w{Ei) has to be uniformly distributed. 

In the case of a Tsallis composition rule one obtains 
E2 — {H - Ei)/{1 + aEi) and we arrive at 



2 woEiiH_Ei} 
d w = -rr—-. — \ ^ dEic 



P {l + aEiY 
This expression can be integrated giving 



d^w = — dpd<p, 
2tt 



(99) 



(100) 



(92) with 



(2 + aH) ln(l + aE,) - ^ {2 + aH + aE,) 
{2 + aH)\n{l + aH) - 2aH 

(101) 

when properly normalized. The only problem is that p{Ei) 
cannot be inverted analytically. Even in the traditional 
case with a = 0, the inversion requires the solution of a 
third order equation: 



p{E, a = 0) = 3{E/Hf - 2{E/Hf 



(102) 



is distributed uniformly between zero and one. It means 
that E is between zero and 77, the total composed energy. 
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Random Energy Density Distribution 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 



Fig. 3. The scaled differential probability density, dw/dx is 
shown as a function of the random energy share of one of the 
collided particles x = Ei/H (cf. ea. (l99p ) for different values of 
aH. Full circles and boxes denote the random energy deviates 
obtained numerically using the rejection method. 



After having Ei and E2 the momenta can be recon- 
structed with the help of a vector triad describing the 
direction of the momentum sum, P. They are given as 
n-i = Pi/P and e; = (— nin2/n±, nin^/n±) where 
the notation n±_ = \Jn\-\- stands for the component 
perpendicular to the first axis. The third orthogonal unit 
vector is fi = (0, —n:i/n±, n2/nj_) . The momentum differ- 
ence vector is hence reconstructed as 



q — Qll n + qj_ (cos (fi e + sin /) 



(103) 



with 



El -El 
2P ' 
Ef + El P^ 
2 4"' 



fl2 



(104) 



3.3 Parton cascade simulation 



First we show some snapshots of the colliding partons in 
the px — Py phase space cut at different stages of the evo- 
lution marked by the average number of collisions per 
particle, t (cf. FigH]). At the beginning i = we pre- 
pared two distributions at a given energy per particle and 
then Lorentz boosted each with ys — 2 units of rapid- 
ity in opposite ways in the Pa;-direction. The dark dots 
represent particle momenta stemming from the respec- 
tively boosted original sets. The evolution towards a zero 
centered and isotropic distribution of momenta signals al- 
ready that thermal equilibration happens. 

Fig [5] presents results of a simple test particle simula- 
tion with the rule h{x, y) = x + y + axy with a — (left) 
and a = 2 (right), respectively. We started with a uniform 
energy-shell distribution between zero and £^0 = 1 with 






Fig. 4. Snapshots of phase space cuts in the px —Py plane for 
colliding partons with the deformed energy composition rule 
h(x,y) = X + y ~ 0.2 xy at t = 0, 0.3, 1, and 3 (from top to 
bottom). 
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Evolution to Boltzmann-Gibbs {a=0) 




Fig. 5. Evolution of single particle energy kinetic energy 
distributions for massless particles towards the Boltzmann- 
Gibbs distribution for h{Ei,E2) = Ei + E2 (left part) and 
towards the Tsallis-Pareto one (right part) for h{Ei,E2) ~ 
El + E2 + 2E\E2- The curves are normalized to the same in- 
tegral / E'^f{E)dE. 



a fixed number of particles N — 10^. The one-particle 
energy distribution evolves towards the well-known expo- 
nential curve for a = 0, shown in the left part of FigO 
These snapshots were taken initially and after 0.1, 0.5, 1, 3 
and 10 two-body collisions per particle. Using the pre- 
scription with a = 2, the energy distribution approaches 
a Tsallis-Pareto distribution. 



4 Non-extensive thermal equilibration 

In order to investigate the equilibration of non-extensive 
systems we start with two subsystems, equilibrated sepa- 
rately. In order to prepare these systems the non-extensive 
Boltzmann equation can be solved numerically in a parton 
cascade simulation as described in the previous section. As 
an alternative way we use initial momentum distributions 
prepared by Monte Carlo rejection techniques in the form 
of ea. (|88p . with different energy per particle but a common 
parameter a for the one an the other half of the particles. 
Then random binary collisions between randomly chosen 
pairs of particles are evaluated. By doing so we apply the 
rules 

X (El) + X {E2) = X (Es) + X (Ei) , (105) 

Pi +P2 P3 +P4- (106) 

In each step of the simulation we select two particles to 
collide. Then we find the value for the new momentum 
of the first particle (ps) satisfying the above constraints 
but otherwise random. Then applying eq. 11061 we calcu- 
late the momentum of the second outgoing particle (P4). 
In these particular simulations we use the free dispersion 
relation for massless particles {Ei{pi) = \pi\), since we 
are interested in the extreme relativistic kinematics case. 
A typical simulation includes 10^ — 10^ collisions among 
10^ — 10^ particles. After 3 — 5 collisions per particle on 
the average, the one-particle distribution approaches its 
stationary form sufficiently. 

The following quantities are conserved during the sim- 
ulation: 

N N N 

i=l i—1 i—1 

(107) 

We use the rule h{x^ y) — x+y+axy for the energy compo- 
sition, here a ~ {q—l)/T is the non-extensivity parameter. 
Our model reconstructs the traditional Boltzmann-Gibbs 
thermodynamics in the limit of a = 0. 

As a preparation for the study of non-extensive ther- 
mal equilibration, we perform simulations on two large 
subsystems with particle numbers A^i = N/2 and N2 — 
N/2, total (quasi-)energies X{Ei) and X{E2) and non- 
extensivity parameter a. The unified system is taken as 
an initial state with iV = A^i + iV2 particles. 



It is in order to make some remark on the energy con- 
servation. For h(x, y) = x + y we simulate a closed system 

with elastic collisions: The sum, Etot — J2iLi-^i' does 
not change in any of the binary collisions. This is differ- 
ent by using a non-extensive formula for h(x,y). With 
a constant positive (negative) a, the bare energy sum is 
decreasing (increasing) while approaching the stationary 
distribution, while the sum of the formal logarithms of the 
energy remains constant. Open systems may gain or loose 
energy during their evolution towards a stationary state. 



4.1 Simulation results 

Our results show that the subsystems do equilibrate, they 
tend towards having a common stationary distribution. 

We present examples with different initial conditions. 
We fix the particle numbers for each subsystem, Ni = 
N2 ~ 100 000. The number of collisions in a typical sim- 
ulation is Ncou = 1 000 000, so that Ncoii/{Ni + N2) = b 
collisions happen per particle. This quantity we use as an 
evolution parameter instead of the real time. This way we 
do not have to know differential cross sections; from the 
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Thermalization of Boltzmann-Gibbs spectra 




Fig. 6. Equilibration of two Boltzmann-Gibbs systems {a — 
0, upper figure) plotted on a linear - logarithmic scale and 
equilibration of two Tsallis-type non-extensive systems (a = 
2, lower figure) plotted on a double logarithmic scale. These 
are results for three-dimensional systems with 10 collisions per 
particle on the average. Each subsystem consists of 100.000 
massless particles. 



viewpoint of the fact of equilibration its rapidity does not 
matter. 

In the figure [6] energy distribution curves are shown: 
the initial and the final ones and the ones after 5 collisions 
per particle, respectively. The upper part plots a Boltz- 
mann system (simulation with a = 0) the lower one a Tsal- 
lis system with the energy composition rule using a = 2. 
In the upper half a logarithmic ~ linear plot is shown while 
in the lower half a double logarithmic plot. These choices 
are selected by the respective high energy asymptotics; 
exponential for a Boltzmann-Gibbs, while power-law for 
a Tsallis-Pareto distribution. It is hard to distinguish the 
energy distributions in the subsystems in the final state, 
the simulation curves are very close to each other. There- 
fore we conclude, that within numerical uncertainties a 
common stationary energy distribution is achieved. 



4.2 Equilibration of large subsystems 

Seeking for a canonical equilibrium state we have to max- 
imize the total entropy given by a general composition 
rule, S{Ei, E2), at the same time satisfying a constraint 
which is in the general case also non-additive: h{Ei,E2) 
is constant. For the moment we neglect the dependence 
on further thermodynamical variables; usually the parti- 
cle number N and the volume V is regarded to be propor- 
tional and extensive. 

In the traditional case both the entropy and the en- 
ergy are combined additively: S{Ei,E2) = S{Ei) + S{E2) 
and h{Ei,E2) = Ei + E2- In the general case by using 
corresponding formal logarithms the quantities Y(S) and 
X{E) have to be considered as additive. Since for asso- 
ciative rules the formal logarithm is strict monotonic, the 
maximum of the total entropy is achieved where Y{S) has 
its extreme. The general canonical principle is therefore 
given by 

Y{S) ~ PX(E) = max. (108) 

The parameter /? at this point is a Lagrange multiplier. 
Applying this for the equilibration of two large subsys- 
tems, and assuming that the entropy of each systems de- 
pends only on its own energy, one arrives at the equilib- 
rium condition 

Y'{S{E,)) Y'{S{E2)) 1 

x'(E,) ^ ^^'^ - x'{E2) ^ ^^'^ r ^^^^^ 

Comparing this with the general canonical form eq. JlOSp 
we obtain that /? = 1/T, and T is an absolute tempera- 
ture in the classical thermodynamical sense. Its relation 
to the entropy, however, has been generalized. In par- 
ticular for an additive entropy, but non-additive energy 
composition rule, one arrives at 1/T = S'{E)/X'{E). 
The relation of this quantity to the logarithmic spectral 
slope, 1/Tsiopc = —dlnf/dE = S'{E) leads to a prac- 
tical tool for the analysis of particle spectra in experi- 
ments. For the Pareto-Tsallis distribution it is given by 
Tsiopc = T/X'{E) = T{l + aE) = T+{q-l)E. The naive 
effort to extract a temperature from energy spectra of par- 
ticles, as it is a widespread usage in relativistic heavy ion 
studies, only works if q = 1, i.e. for spectra exponential in 
the particle energy. Otherwise an energy dependent slope, 
and a curved spectrum in the logarithmic plot has to be 
interpreted. 

The inverse logarithmic slopes of single-particle kinetic 
energy spectra in the generalized case are functions of the 
energy: 

Tslopc = , . ■ (110) 

For the Tsallis-Pareto distribution they are linear func- 
tions, Tgiope = T + (q — 1)E. Such slope parameters are 
plotted in Ref.[55] for the respective subsystems before 
and after equilibration (10 collisions per particle on the 
average). Within numerical uncertainties it is clear that 
common-a systems do equilibrate at a common tempera- 
ture also in the a ^ case. 
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Fig. 7. The evolution of the Boltzmann entropy per particle 
during collisions with non-additive energy composition rules: 
the hotter body cools, the cooler body warms up, while the 
total entropy also increases. In the insertion a magnification of 
the curves is shown. 

The (in our case Boltzmann) entropy also evolves due 
to the collisions. In Fig[7]the evolution of the entropy per 
particle is plotted for the hot and cool subsystems, and 
for the total system respectively. Since the composite sys- 
tem is combined from equal numbers of particles in each 
subsystem, the total entropy per particle starts with the 
arithmetic mean of the respective specific entropies. This 
value, however, rises somewhat, featuring a trend accord- 
ing to the second law of thermodynamics. 

5 Power-law tailed hadron spectra, flow and 
quark coalescence 

As an application of the above reviewed treatment of non- 
extensivity, in this section we demonstrate that hadronic 
transverse momentum spectra stemming from relativistic 
heavy ion collisions can be well described by cut power- 
law spectra in statistical models. In order to do so, one 
has to disentangle effects of a possible transverse flow on 
these spectra and then test whether the result complies 
with the thermal assumption; i.e. that the dependence on 
momenta is through a dependence on the kinetic energy, 
E ^ fi = E — m only. The fi ~ m assumption corresponds 
to a vanishing Fermi momentum for fermions, so this is 
the natural assumption at zero net baryon density. There- 
fore transverse momentum spectra at mid-rapidity are ex- 
pected to follow such statistical model assumptions the 
best. 



5.1 Spectral temperatures in relativistic heavy ion 
collisions 

It has been long discussed, how a temperature can be con- 
jectured from observations on particle spectra produced 
in relativistic heavy ion collisions. One intriguing way is 



to look at the transverse momentum, px, spectra around 
mid-rapidity. The different identified hadrons, mostly pi- 
ous, kaons, protons and antiprotons, have to demonstrate 
that their abundance in the momentum space depends on 
their kinetic energy; this phenomenon at zero rapidity is 
the so-called — m-scaling. The transverse mass is given 
as niT = \/ rn^ + p^j strictly zero rapidity this is the 
total relativistic energy. 

The analysis is made a little more involved by the fact 
that the source emitting the detected hadrons is flowing 
in all directions. The most prominent effects are due to a 
relativistic transverse flow with velocity vt (and a corre- 
sponding Lorentz factor 77^ = 1/ ^yl — v'^ in units where 
c = 1). The relativistic energy of a particle in the frame of 
the emitting source cell is given by the Jiittner variable: 

E ~ u^p^ — ^tttit cosh(y — 77) — jt^tPt cos{(p — <P). 

(Ill) 

Here the four-velocity of the source and the actual four- 
momentum of the particle are parametrized by rapidity 
and angle variables: 

Up. = (7t coshr;, 7t sinhyy, 7tZ't cos<?, 7tWt sin<?), 
— {rriT cosh y^rriT sinh y,pT cos if, pt sill ip). (112) 

We consider a thermal model for the particle spectra; then 
the yield is supposed to depend on the Jiittner variable E 
given by eq. ljllip . Assuming a general distribution f{E) ~ 
exp{—{X(E) —m)/T), which is monotonic decreasing, one 
finds its maximum at the minimum of E. This variable is 
minimal at the rapidity ymm — ^-i^d angle t/Smin = 
giving 

Emin = irniT - ItVtPt- (113) 

This Lorentz-boosted transverse energy reaches its mini- 
mum at the transverse momentum value pT.min = rwyTVT, 
leading to rriT.min = m'^T and -Emin — m- The expansion 
around this minimum in the px-distribution is an effective 
Gaussian: 

e-i--y-^e.J-^^^^^). (114) 
\ ZrwyT -f 7T / 

In fact, according to experimental findings at RHIC 
the observed particle spectra have to be corrected for a 
transverse flow in order to reach mr-scaling. 

5.2 Non-extensivity in quark matter and in hadron 
matter 

We conjecture that the power-law tails observed in hadronic 
spectra may stem from non-extensivity of the suddenly 
hadronizing quark matter. We look for a connection be- 
tween quark and hadron spectra in the framework of the 
quark coalescence model. A coalescence of a quark and an 
antiquark into a meson produces a yield proportional to 
the quantity: 

F{p)^ J f{E{P/2 + q))f{E{P/2~q))C{q)d^q. 

(115) 
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HADRON FITTINGS 




CO-MOVINO MOMtN I LIM 




F(P)« r{E{P/2)) I C{q)dPq. 



(116) 



In particular light hadrons made from massless quarks fol- 
low the quark-scaling rule: 



h^droniE) oc r{E/n). 



(117) 



As a consequence particular properties of the non-extensive 
thermal model between quark and hadron matter also 



scale: T, 



mesons 



Tbaryons = ^quarks for the temperature, 

while 9mesons-l = (^quarks " l)/2 for meSOUS and gbaryons- 

1 — (?quarks — l)/3 for baryons. Since for a Tsallis-Pareto 
distribution the inverse logarithmic slope turns out to be 



slope 



^T+iq-l)iE~m), 



(118) 



the rise of these slopes reflect the non-extensivity param- 
eters. In FigElwe show the test of the coalescence model 
prediction for the meson to baryon ratio. 




Fig. 9. Inverse logarithmic slopes, Tsiopc(£) ~ —dE/d\nf{E) 
extracted by numerical derivation from the experimental 
hadronic spectra (after subtracting a common flow effect). The 
full lines correspond to a common meson and baryon fit their 
steepness keeping the ratio 2:3 predicted by the quark coales- 
cence picture. 

Extracted values for the quarlt q-parameter 



Fig. 8. General shape of pr spectra for pions, kaons and an- 
tiprotons in relativistic heavy ion experiments (upper figure). 
A fit is done by using for X{E) the Tsallis-Pareto form with 
parameters T and a, corresponding to a common temperature 
of T{mi) = 0.160 MeV for the different particles, and a trans- 
verse flow velocity vt = 0.52. In the lower part the ratio of 
the Tsallis fit to the experimental values can be inspected in a 
linear plot. 



Here we integrate over the relative momentum of the quarks 
with a coalescence factor, C(q), for which a simple model 
has been utilized |23j . For common momenta much larger 
than the relative one \P\ ^ |q| on obtains 




Quark coelasecence slopes at zero p, 





0.21'm+0.052 
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Fig. 10. The q parameter of quark matter extracted from 
hadronic spectra assuming quark coalescence at a sudden 
hadron formation (upper picture). The spectral inverse slope 
as a function of the minimal = m agree with the 

linear prediction from the coalescence scaling. 



Furthermore these predictions of the non-extensive phe- 
nomenology meet the curves from pQCD calculations, with 
the following surmised properties of quark matter at RHIC: 
T = 140 . . . 180 MeV, q = 1.22, vt = 0.6. [Mj- Deviations 
from the product rule suggested by the simplest quark co- 
alescence idea occur at pr values lower than 1 GeV. We 
attribute these to a further constituent in real hadrons, 
namely non-perturbative gluons simulated by a string en- 
ergy contribution 26] . 

We note that a stringy interaction remainder above 
the color deconfinement temperature Tc ~ 170 MeV in 
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quark gluon plasma also describes the main effects on the 
quark matter equation of state seen in lattice QCD calcu- 
lations successfully [27J. Both the presence of a string like 
pair potential for a however minor percentage of pairs, as a 
microscopic model, and the assumption of non-extensivity, 
g > 1 , as a descriptive phenomenology are able to explain 
the value of the energy per particle, E/N = 6T = 1 GeV, 
which has been found by fitting statistical hadronic reso- 
nance gas models with Boltzmann distributions. 

Finally some remarks are in order to the Tsallis-Pareto 
fits to energy spectra. In several cases naively a fit is done 
in the original form: 

/(i^)- (l + (9-1)1) (119) 

to which the following inverse logarithmic slope depen- 
dence belongs: 

T.10PC = T + (g - 1)£;. (120) 

A more sophisticated approach (as one suggested in pT]). 
however, uses the original Tsallis-Pareto form for the num- 
ber density distributions of particles and for the generat- 
ing thermodynamical potential, for the logarithm of the 
canonical partition function. This way in this second ap- 
proach the energy distribution is described by the q-ih 
power of the naive factor: 

/(i^)- (l + (9-l)f) (121) 
and the corresponding inverse slope 

T.iopo = -.f+(^--}j E. (122) 

We observe that the qualitative behavior is the same, but 
the interpretation of the fit parameter is different in these 
different approaches. The correspondence between the en- 
ergy spectrum fit parameters is given as 

T = f/q, 

q = 2- 1/q. (123) 

In a sense q and q are double-duals of each other, both 
using the 1/q- and the 2 — g-duality. Also the estimated 
temperature parameter differ. Typical values from rela- 
tivistic heavy ion experiments are q ^ 1.2 and q « 1.25, 
as well as, f = 1.25T. 



6 Conclusion 

In conclusion we reviewed basic concepts of non-extensive 
thermodynamics which may be relevant in understanding 
particular features of hadronic spectra stemming from rel- 
ativistic heavy ion collisions. The overall presence of rel- 
ativistic speeds of particles in the physical system under 
investigation on the other hand offers a unique possibility 



to study and - whenever necessary - to generalize familiar 
thermodynamics. 

We presented some general arguments for a possible 
need to face with total energy and entropy not being 
proportional to the particle number even in the large N 
limit. These arguments are based on the long range na- 
ture of pair interactions. This phenomenon, called non- 
extensivity, was then related to the generalization of com- 
position rules of the familiar thermodynamical extensives, 
like energy and entropy. We have mathematically proved 
that abstract composition rules become symmetric and 
associative in the large N limit, provided that the compo- 
sition function, h(x, y), is at least right-sided differentiable 
at y = 0"*". This means that associative composition rules 
constitute attractors among all rules when approaching 
the thermodynamical limit. As a consequence the associa- 
tivity of the composition rule is a thermodynamical re- 
quirement. 

The key quantity in this proof, the formal logarithm, 
relates the abstract composition rule to the addition of 
the system size indicator, to the particle number N. We 
gave the formula how to construct it. Based on the formal 
logarithm the widely used deformed exponential and log- 
arithm functions can easily be derived. While the former 
describes the energy distribution in canonical equilibrium, 
the latter defines a generalized formula for the entropy. An 
additive entropy can be always gained from this expression 
by taking its formal logarithm. 

We presented some often used composition rules to- 
gether with the corresponding equilibrium energy distri- 
butions and entropy formulas including the traditional 
Boltzmann-Gibbs formula derived from the simple addi- 
tion (extensivity) , the Tsallis rule, the Kaniadakis rule 
and - for the sake of demonstration - the Einstein rule 
for composing relativistic velocities. Our general method 
in this case leads to the rapidity as the additive formal 
logarithm. Among non-associative composition rules the 
class of h{x,y) = x-\-y-\-G{xy) is found to be particularly 
interesting in high energy physics, since it asymptotically 
approaches the Tsallis rule leading to power-law tailed 
energy distributions in canonical equilibrium. We have 
demonstrated that such a composition rule may emerge 
in the extreme relativistic kinematics limit from an en- 
ergy correction to a pair of particles in a medium which 
is a function of the Lorentz invariant relative momentum 
squared variable Q^. In fact this is frequently the case 
when following several interactions among partons accord- 
ing to the formulas derived from (or at least motivated by) 
QCD. Finally it is interesting to note that the elementary 
property, h{x,Q) — x is related to (t(1) = property of the 
entropy density function if the composition rule h{x,y) 
is assumed for composite states with factorizing proba- 
bilities. In this case the general result a{p) = lna(l/p) 
emerges in the thermodynamical limit, with In^ = oln 
being the corresponding deformed logarithm function. 

The - in some sense opposite - requirement, i.e. aiming 
at an additive entropy formula while the probabilities do 
not factorize, but their logarithms follow a general com- 
position rule h{x,y) instead of the usual addition, leads 



16 



Tamas S. Biro, Gabor Purcsel, Karoly Urmossy: Non-Extensive Approach to Quark Matter 



to a more complex relation between the formal logarithm 
of the rule, L, and the entropy density a. 

In the second part of this review we compiled the most 
important numerical results on parton cascade simulations 
of non-extensive systems. Following the presentation of 
a class of generalized Boltzmann equations, and proving 
that the second law of thermodynamics can only be ful- 
filled if the derivative of the entropy density is a linear 
expression of the deformed logarithm of the one-particle 
phase space density (cf. eq. (|80|) ). we presented some de- 
tails of the kinematical description of relativistic particles 
in such simulations. We payed special attention to the ran- 
dom choice of particle momenta after an, in energy non- 
additive, pair collision (which can have a physical reason 
in the influence of third or further particles, or fields in a 
dense medium). 

Results on the phase space evolution under non-exten- 
sive energy composition rules were presented demonstrat- 
ing the ability of such a computer simulation to generate 
power-law tailed energy spectra in the detailed balance 
state of the non-extensive Boltzmann equation. The im- 
portant question of equilibration between two large sub- 
systems, related to the zeroth theorem of thermodynam- 
ics, was also investigated by us numerically in this frame- 
work. We found that non-extensive systems with the power- 
law tailed Tsallis-Pareto energy distributions do behave 
as they should, just the thermodynamic temperature, T 
is related to the microcanonical equation of state, S{E), 
by receiving corrections due to the formal logarithms of 
the entropy and energy composition rules (cf. eq. (|109p ). 

Finally our studies on the hadronization of quark mat- 
ter in relativistic heavy ion collisions revealed that if the 
quark coalescence is a dominant mechanism, then the non- 
extensivity parameter, q — 1 = aT also must show the 
quark number scaling. This assumption can be and should 
be tested on experimental data and should be related to 
other information on quark coalescence, e.g. to those ob- 
tained from studies of the elliptic flow. 

Certainly there remain open questions for further re- 
search. Among them the study of the quark matter equa- 
tion of state with elementary field theory means, as lat- 
tice QCD, in a non-extensive canonical state is still a hard 
challenge. Also the determination of the pair correlation 
function, g{r), from first principles in microscopical cal- 
culations should help to identify those physical situations 
where the concepts and formalism of non-extensive ther- 
modynamics have to be used. Meanwhile the physical rea- 
son for a non-exponential energy distribution can be nu- 
merous. The quark gluon plasma is a wonderful candidate 
for finding non-extensive behavior, since long range effects 
are there at any finite temperature. 

Acknov^rledgments 

This work has been supported by the Hungarian Na- 
tional Science Fund, OTKA (K49466, K68108). Discus- 
sions with C. Tsallis, G. Wilk, T. Kodama, G. Kaniadakis, 
P. Van and A. Laszlo are gratefully acknowledged. 



References 

1. Shannon C. E., Bell Syst. Tech. J. 27(1948)379; ibid. 
27(1948)623. 

2. Jaynes E. T., Phys. Rev. 160(1957)171; ibid. 
168(1957)620. 

3. Renyi A., Acta Math. Acad. Sci. Hung. 10(1959)193; 
Probability Theory, North Holland, Amsterdam, 1970. 

4. Daroczy Z. Aczel J., On Measures of Information and 
their Characterization, Academic Press, New York, 1975. 

5. Fisher R. A., Statistical Methods and Scientific Inference, 
Oliver and Boyd, Edinburgh-London, 1959. 

6. Van P., Physica A365(2006)28. 

7. A. Renyi, Probability Theory, North Holland, Amster- 
dam, 1970; A. Wehrl, Rev.Mod.Phys. 50, 221, 1978; 
Z. Daroczy, Inf.Control 16, 36, 1970; J. Aczel, Z. Daroczy, 
On Measures of Information and their Characterization, 
Academic Press, New York, 1975. 

8. C. Tsallis, J. Stat. Phys. 52, 50, 1988; Physica A 221, 
277, 1995; Braz.J.Phys. 29, 1, 1999; P. Prato, C. Tsal- 
lis, Phys.Rev.E 60, 2398, 1999; V. Latora, A. Rapisarda, 
C. Tsallis, Phys.Rev.E 64, 056134, 2001; Physica A 305, 
129, 2002. 

9. Biro T. S. Purcsel G., Phys. Rev. Lett. 95(2005) 162302. 

10. T. S. Biro, Eur. Phys. Lett.xx(2008)pp. e-Print: 
arXiv:0809.4675 [nucl-th] 

11. Tsallis C., Gell-Mann M. Sato Y., Proc. Nat. Acad. 
Sci.lO2(2005)15377. 

12. Castillo E., Iglesias A. Rm'z-Cobo R., Functional Equa- 
tions in Applied Sciences, Elsevier, 2005. 

13. Tsallis C, J. Stat. Phys.52(1988)479. 

14. Kaniadakis G., Phys. Rev. E66(2002)056125; ibid. 
72(2005)036108; Physica A365(2006)17. 

15. G. Kaniadakis, Physica A 296, 405, 2001; Phys.Rev.E 66, 
056125, 1, 2002. 

16. J. A. S. Lima, R. Silva, A. R. Plastino, Phys.Rev.Lett. 
86, 2938, 2001 

17. A. Lavagno, Phys.Lett.A 301, 13, 2002. Relativistic 
Nonextensive Thermodynamics, 

18. C. Tsallis, E. P. Borges, cond-mat/0301521| C. TsaUis, 
E. Brigati, cond-mat/0305606 ; C. Tsallis, Braz.J.Phys. 
29, 1, 1999; A. Plastino, A. R. Plastino, Braz.J.Phys. 29, 
50, 1999. 

19. Q. A. Wang, A. Le Mehaute, J.Math.Phys. 43, 5079, 
2002; Q. A. Wang, Chaos, Solitons and Fractals 14, 765, 
2002; Eur.Phys.J.B 26, 357, 2002. 

20. S.Abe, Physica A 300, 417, 2001; Phys.Rev.E 63, 061105, 
2001. 

21. L. Borland, Phys.Rev.E 57, 6634, 1998; D. H. Zanette, 
Braz.J.Phys. 29, 108, 1999; G. Kaniadikis, Phys.Lett.A 
283, 288, 2001. 

22. Biro T. S. Purcsel G., Phys. Lett. A372(2008)1174. 

23. Biro T. S., Levai P. Zimanyi J., J. Phys. G28(2002)1561, 

24. Biro T. S., Urmossy K. Barnafoldi G. G. , J. Phys. 
G35(2008)044012, 

25. Biro T. S. Urmossy K, Eur. Phys. J. ST155(2008)1. 

26. Biro T. S. Urmossy K, Contribution to SQM2008, held 
at Tsinghua University, Beijing, China, Oct. 6-10. 2008 

27. Biro T. S. Cleymans J., Phys. Rev. C78(2008)034902. 



